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The perturbative QCD formula for minijet production consitutes an important ingredient in 
models describing the total cross section and multiparticle production in hadron-hadron scattering 
at high energies. Using arguments based on s-channel unitarity we set bounds on the minimum 
value of pt for which the leading twist minijet formula can be used. For large impact parameters 
where correlations between partons appear to be small we find that the minimum value of pj should 
be greater than 2.5 GeV for LHC energies and greater than 3.5 GeV for cosmic ray energies of 
about 50TeV. We also argue that for collisions with values of impact parameters typical for heavy 
particle production the values of minimum pt are likely to be considerably larger. We also analyze 
and quantify the potential role of saturation effects in the gluon density. We find that although 
saturation effects alone are not sufficient to restore unitarity, they are likely to play an important 
role at LHC energies. 
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I. INTRODUCTION 

The high energy hadronic interactions which will be 
studied at the LHC, and the interaction of ultra-high 
energy cosmic rays in the atmosphere both involve 
QCD effects that are still poorly understood and largely 
unexplored experimentally. Therefore, in extrapolat- 
ing Monte Carlo simulations to these very high ener- 
gies, one often must resort to phcnomenological mod- 
els p], 2, 3, 4, 5, 6] based on fits to lower energy data. Of 
particular interest is the contribution from jets with mod- 
erate relative transverse momentum pt (so-called mini- 
jets or semi-hard jets). 

The basic input to these models is the total inclu- 
sive minijet cross section in hadron-hadron scattering, 
typically calculated using the leading-twist perturbative 
QCD (pQCD) factorization formula. Schematically, the 
minijet cross section is the convolution integral of a par- 
ton density for each colliding hadron and a partonic cross 
section for two-to-two parton scattering, 
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fh2{x2,^J-'^), (1) 



where (8) denotes convolution integrals in longitudinal 
momentum fractions Xi,X2 for each of the colliding 
hadrons. The scale fi is the factorization scale which is 
usually set to be equal to pt, the transverse momentum. 
Equation ([T]) has to be regulated in the low momentum 
regime by introducing a cutoff p^. Equation ([1]) consti- 
tutes the semi-hard part in most models of high energy 



hadron scattering. The part which includes soft interac- 
tions must be modeled separately. 

Since the integrand in Eq. ([1]) is sharply peaked at 
very low transverse momenta, the minijet cross section is 
highly sensitive to the value of the cutoff p^. The problem 
of determining a minimum pt is generic to all the models 
which aim to describe the particle production in hadron 
collisions Ili,ld,l3|,y,l5|,l6!,l3,|8|. 

The true range of validity for pi remains a point of 
some controversy. The leading twist pQCD expression 
for Eq. ^ is most reliable when p^ is large. However, 
one hopes to utilize the full power of pQCD for as wide 
a range of kinematics as possible. 

Imposing unitarity in impact parameter space provides 
a potential constraint on allowed values of p^ at very 
high energies. Unlike in the case of the total cross sec- 
tion, however, the unitarity constraints cannot be im- 
posed directly on ([1]), which is an inclusive cross section 
proportional to average jet multiplicity. Therefore, less 
direct methods which take jet multiplicity into account 
are needed. A parton model based picture of multiple in- 
teractions has been shown to be a natural framework Q 
with which to address the unitarity issue. This approach 
has been applied, for example, in (lo| in the context of 
Monte Carlos for cosmic ray air showers. 

In this paper, we propose a simple method for con- 
straining the minimum value of p^ at large impact pa- 
rameters (6 > 1.5 fm) for which the minijet formula in 
Eq. Il]) makes sense. Our method is based on s-channcl 
unitarity of the amplitude in the impact parameter repre- 
sentation. Our results are sensitive only to the transverse 
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distribution of gluons in the nucleon which is measured 
in small- a; hard exclusive processes. We produce conser- 
vative bounds on - the true minimum should be 
higher than what we find. In this paper we restrict at- 
tention to large impact parameters where it is reasonable 
to neglect parton correlations. Apart from this, however, 
we avoid discussing any specific model of the underlying 
dynamics responsible for unitarization. Nevertheless, the 
approach we use in this paper is likely to be applicable 
at smaller impact parameters, given a particular model 
of parton correlations. 

The lower bound we obtain on the transverse momen- 
tum cutoff is at least 2.5 GeV for center-of-mass pp colli- 
sions at 14 TeV and impact parameters > 1.5 fm. We also 
consider here much higher energies, -^/s « 50 TeV which 
are relevant to cosmic ray studies. Here, we find that the 
minimum cutoff probably becomes closer to 3.5 GeV for 
the same values of impact parameters for pp collisions. 
We also discuss possibilities for extending the analysis to 
smaller impact parameters. 

The structure of the paper is as follows: in the next 
section we introduce the formula for minijet production. 
By investigating the integrand at fixed values of pt and 
the energy ^/s we constrain the range in the rapidities 
and X values that dominates the integral. On the ba- 
sis of the saturation scale from the Golec-Biernat and 
Wusthoff model (GBW) we estimate the role of the sat- 
uration corrections in the gluon density at the Tevatron 
and at the LHC. 

In Sec. Ill we develop a general method based on uni- 
tarity arguments for constraining the minimum value of 
the Pt for the minijet formula to be used. 

In Sec. IV we give a general discussion of the radius of 
the interaction in the context of the eikonal model. We 
make comparisons here to the interaction radius for deep 
inelastic scattering (DIS) 

Finally, in the appendix we include the formulae for the 
dependence of the mass parameter on x and the scale in 
the profile function. 

II. MINIJET CROSS SECTION 

We begin by writing Eq. ([T]) more explicitly. The inclu- 
sive minijet cross section in hadron-hadron scattering is 
typically calculated using the leading-twist perturbative 
QCD (pQCD) factorization formula, 

"■Pt 

Here, the differential cross section da /dpi is for 
two-to-two parton scattering. In this work we calculate 
it at the lowest order in pQCD. The standard leading 
twist parton distribution function, fj{x,^'^) is evaluated 
at a hard scale /i^ which is typically chosen to be the 



transverse momentum squared, p^. The factor 1/(1 -I- ^fc;) 
takes into account necessary symmetry factors when the 
outgoing partons are identical. The if is a factor used to 
correct for higher order terms, and is fixed to 2 in many 
models. In NLO treatments of single jets, K is found to 
be closer to 1.0. However since processes where at least 
three jets are produced appear to be dominant at these 
energies K — 1 corresponds to effective K 1.5 for 
our observable. In keeping with our goal of maintaining 
conservative estimates of unitarity bounds, we take this 
value for the K factor. 

The kinematic limit for the xi,X2 integrations is 

XIX2S > A{p^f . (3) 

As we will see shortly, in the UHE limit and at small pt , 
the dominant range of the integration in ([2]) is over small 
values of x where the gluon density, g{x,p1) dominates. 
A continuing source of uncertainty is the lowest value of 
the for which Eq. ([2]) makes sense. Since the cutoff is 
in the regime where the integrand begins to spike, small 
changes in the cutoff can lead to dramatic changes in 
the integrated cross section, Eq. ([2]). This is illustrated 
in Fig. [T] where the integrated minijet cross section is 
shown as a function of the c.m.s. energy for different 
values of the p^. We see a drastic change in the nor- 
malization, although the shape and the rise of the cross 
section with the energy is universal. This sensitivity to 
is a symptom of our limited understanding of physics in 
transitional region between small values of pt, where the 
distribution in transverse momentum must be modeled, 
and large values where Eq. ([2]) is a good approximation. 

A. The Kinematic Range of Validity for the 
Minijet Formula 

In order to gain an intuition about the structure of the 
minijet cross-section at very high energies, it is useful 
to study the integrand of Eq. ([2]) as a function of the 
rapidities 2/1,2/2 of the jets. We show a contour plot in 
Fig. [T] for the integrand versus j/i and y2^ related to xi 
and X2 via the equation. 

We use CTEQ6M [l^ parton distribution sets for eval- 
uation of Eq. The smallest contours represent the 
maxima of the integrand. To get an idea of the important 
regions of the integrand, we have drawn solid thick-red 
curves along the contours marking half of the peak value. 
Plots (a) and (b) show the result of calculations done at 
1.96 TeV and 14 TeV respectively. In both cases, pt is 
fixed at 2.5 GeV. Note that as ^/s increases, the peaks 
begin to form long ridges that extend toward the lower- 
left and upper-right corners of the plots where yi and 2/2 
have the same sign. The ridges also expand if we decrease 
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Pt, as shown in plot (c) where pt = 1.5 GeV, still with 
= 14 TeV. 

For figure m^a) the peaks are located at rapidities, 
2/1 ^ -y2 ^ ±1-2, 

and so the x values probed in the collision with this con- 
figuration are approximately equal a;i ~ a;2 ~ 0.005. 
Still, at rapidities 

lyil - 4, 12/2I ^ 2, and signyi = signj/2 , 

the integrand is rather large. Such a configuration corre- 
sponds to asymmetric values of xi,X2, 

xi,2 =i 2 X 10""* < X24 ~ 0.08 . 

At the maximum LHC energy - 14 TeV - plotted in figure 
(b), the central values are, =4> Xi <-*■ X2 

xi~X2^ 0.0008, 

but the ridges extend toward values of rapidity, 

|yi|~6,|2/2|~4, 

corresponding to xi^2 — 4 x 10~^ <^ X2,i — 0.08 . There- 
fore, both aligned and anti-aligned jets become impor- 
tant. 

In the regions where yi and 7/2 are both very large 
and of the same sign, one of the values of x becomes 
very small. The is the regime where one of the giuon 
densities begins to blow up, and one expects non-linear 
effects (saturation) to come into play, see for example 
[lll)[i2l- For a given value of x saturation models pre- 
dict that g(x, p) becomes saturated for values of /i below 
the saturation scale, Qs- Studies of deep inelastic scat- 
terin g using various versions of the pQCD dipole pic- 
ture [il, El EKES have resulted in rough bounds on 
where saturation effects should begin to become signifi- 
cant. To illustrate the sensitivity of the integrand to the 
saturation effects in the gluon density we superimpose the 
saturation scale from the GBW model [l3|,[l3| onto our 
contour plots. The blue lines in the corners of the con- 
tour figures correspond to the saturation line [2^ from the 
GBW model. These lines mark off the regions of satura- 
tion in the upper right-hand and lower left-hand corners 
of the figures. Values of 2/2 that lie within the blocks 
marked off by the saturation line indicate regions where 
either g{xi,pt) or g{x2,Pt) is likely to be very sensitive to 
saturation effects since either Qs{xi) or Qs{x2) are larger 
than the fixed value of pt- The percentage of the integral 
affected by the saturation corrections is about ~ 11% 
for the Pt = 2.5 GeV and - 53% for pt = 1.5 GeV for the 
LHC energy. We did not show the saturation lines for the 
Tevatron energy since for this choice of pt — 2.5 GeV the 
saturation effects are very small. We note, however, that 
the model used in [l^ allows xi,X2 to go to somewhat 
smaller values without actually violating unitarity. 

When evaluating the saturation effects one has to take 
into account the difference between the saturation scale 
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FIG. 1: Inclusive minijet cross section evaluated from ([2} for 
three different cutoff values pi = 1.5 (dotted line), pi = 2.5 
(dashed line) and pi — 3.5GeV (solid line). The K factor is 
set to one in the plots. Parton distributions are CTEQ6M. 

for the quark and gluon dipolcs. Since here the domi- 
nant process is gg — s- gg scattering the saturation scale 
is enhanced by an additional color factor Ca/Cf = 9/4 
relative [l^ to the case for the quark scattering. We 
therefore include the corresponding rescaled saturation 
region in Fig. [21 marked off by the solid blue line (as 
opposed to the dashed-blue which corresponds to the 
original unrescaled value of the saturation scale from the 
GBW model). 

Furthermore, our analysis has only considered the lead- 
ing gg —^ gg contribution so far. In the calculation of 
the total inclusive minijet cross section we include con- 
tributions from all channels. The gg —^ gg channel is of 
course dominant, it is about 60 — 70% of the total value 
depending on the choice of the parton density and the 
c.m.s. energy. The qg — s- qg contribution however, is 
non-negligible and constitutes about 30% of the minijet 
cross section. It is particularly important for asymmetric 
configurations of the minijets. In this case the large con- 
tribution comes from the region where the quark density 
is evaluated at large x and the gluon density is evaluated 
at very small values of x. 

Figures [2{a-b) demonstrate that gluon saturation ef- 
fects will likely be significant at LHC energies and higher. 
However, we stress that other effects (e.g. multiple scat- 
tering) are equally important for unitarizing the pp cross 
section. In fact the multiple hadron scattering effects 
are likely to become dominant at the lower energies, in a 
regime where the saturation of the gluon density is not 
very large. In this paper, we take a very general ap- 
proach that does not distinguish the underlying physics 
responsible for unitarization. 



III. IMPACT PARAMETER DEPENDENCE OF 
THE MINIJET CROSS SECTION 

In this section we review the basic elements needed to 
discuss unitarity in impact parameter space. 
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(c) 

FIG. 2: Rapidity distribution of the integrand in the minijet cross section formula for (a)^^ = 1-96 TeV, pt ~ 2.5 GeV, (b) 
-^5=14 TeV, Pt = 2.5 GeV and (c) ^/s=14 TeV, pt = 1.5 GeV. The calculation is done for the gluon-gluon channel only. The 
thick - red contour is where the integrand is 50% of the maximum value. The corner regions marked off by the dotted blue lines 
indicate where saturation effects may be present according to the GBW model. The solid blue line is the result of rescaling the 
saturation scale by a factor of 9/4 as discussed in the text. 



A. Basic Formulation 

Following standard steps the total, elastic and inelastic 
cross sections can be written at high energies as, 

(Jtotis) = 2 J d^hReris,b), (4) 

<Jei{s) = J d^h \T{s,b)\\ (5) 

a^nelis) = ^ d^b (2 Rc r(s, 6) - |r(s, 6) | , (6) 

where the profile function T{s,b) is the Fourier transform 
in impact parameter of the scattering amplitude for two- 
to-two process 

r(M) = ^-l^/d^qe--A(M). (7) 



Assuming that the amplitude is dominantly imaginary 
(which is appropriate at high energies), the unitarity con- 
straint on r(6, s) is, 

r(s,6)<i. 

B. Impact Parameter Dependence for Hard 
Scattering 

We appeal directly to experimental data to obtain the 
impact parameter dependence of the hard collisions. We 
use the generalized gluon distribution function extracted 
directly from J/^E* electroproduction production [20l |. 

The generalized gluon distribution is related to the 
standard gluon distribution function and a gluon form 
factor Fg(x, t, jj?') through the following defining relation 
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where 



Fgix,t = 0,^^) = 1 . 



(9) 



The Fourier transform of the ghionic form factor gives 
the profile in impact parameter space 

(10) 

where the integration over p is 2-dimensional. We have 
the normahzation condition 



(PpTg{x,p,p) = 1 , 

which is a trivial consequence of the previous condition 
([9|). Interpreting J-'g{x, p, fi) as the transverse spatial 
spread of hard partons, we may write Eq. Q in the form, 

where we have defined, 

J^2{b,s,pl) T ^ [ dxidx2 [ dp^x 

^ 1 + Sk.i Jo J 

fg{x2,pf)P2{b,Xi,X2,Pt)0{pt -pD, (11) 
with an over-lap function given by, 

P2{b,Xi,X2,fJ.) = 

J d'^piJ^g{xi,\pi\,p)J^g{x2,\h-px\,p). (12) 

As usual, we use p, ~ pt for the hard scale. It follows from 
these definitions that A/'2(fo, s,pj) should be interpreted as 
the relative probability for producing at least one minijet 
pair at impact parameter b with pt > p^ and center of 
mass energy ^/s. 

For a wide range of energies, a reasonable parameteri- 
zation of the two-gluon form factor is the model of Frank- 
furt, Strikman, and Weiss (FSW) [U, 



Fg{x,t,p) 



1 



(13) 



Spreading of the form factor at small- a; and evolution 
with the hard scale p are taken into account by allowing 
mg to vary with x and /z. For this, we use the param- 
eterization given in [2ll | (see Appendix B). Performing 
the Fourier transform in Eq. (jTU]) while using yields 
an analytic expression for the impact parameter space 
gluonic form factor, 



Tg{x,p,p) 



ml{x,p)p 



We make one further approximation to simplify the anal- 
ysis. The dependence of mg on x is rather slow. There- 
fore, we make the replacement, 

P2{b,Xi,X2,Pt) > P2{b,x,x,p1) = P2{b,s,p1), (15) 

where x = 2p'i/^. We have calculated P2{b,xi,X2,Pt) 
numerically, and have verified that it differs from 
P2{b,s,Pf) by no more than ten percent within the es- 
sential region of the integrand in Eq. (jlip. We can then 
approximate, 

AA2(6, « a%l,{s,p'i)P2{b, s,p'i) . (16) 

P2(&, s,P() can then be determined analytically from 
Eq. nil): 

"^g(^,P?) ( mg{x,p'i)b^ ^ 
127r 



P2ib,s,p't) 



K3{mg{x,p^t)b). 

(17) 

This quantity is normalized to unity Jd^bP2{b,s) = 1 
and represents the probability density in impact param- 
eter space for the minijet production. Thus, the average 
squared impact parameter for minijet production is. 



'b') = 



J d'h b'P2{b,S,p'i). 



(18) 



This provides a quantitative measure of the width of the 
overlap function which can then be used to compare with 
other models. 



IV. UNITARITY CONSTRAINTS 

We address the unitarity issue directly by maintaining 
the impact parameter representation and by calculating, 
within a particular model, the relative contribution to the 
total inelastic cross section from the exclusive dijet and 
2fc-jet production cross sections. In this section, we start 
with a simple model with uncorrelated jet production at 
large impact parameters. 



A. Multiple Collisions 

We follow steps similar to Ref. Q to take into account 
multiple collisions. In the simple parton model picture, 
each collision is a two-to-two collision between partons 
resulting in an outgoing minijet pair. For production of 
k minijet pairs, with fc > 1, there is a convolution integral 
analogous to Eq. (fTT|) which we symbolize by N2k {b, s) . 
(Multiple jet production will also involve cutoffs on trans- 
verse momentum.) For k minijets, the normalization con- 
dition is, 



d^bU2k{b,s) 



<^2k [S) 



An 



Ki{mg{x,p)p). (14) 



where a^'^ is the 2A;-jet inclusive cross section. We will 
refer the Af2kib,s) as "impact parameter cross sections" 
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since they integrate to total inclusive cross sections for 
producing 2k minijets. We have not used the symbol 
r(s, b) since this is already used to represent the profile 
function Eq. ^ which integrates to the total cross sec- 
tion. Next, we define the Af2k{b,s) {k > 1) to be the 
exclusive analogues of the N'2k{b, s). That is, 



is the cross section for producing exactly k minijet pairs 
at impact parameter b. From the definition of exclusive 
quantities, the inclusive 2-jet impact parameter cross sec- 
tion can be expressed in terms of the exclusive impact 
parameter cross sections by writing, 



Af2{b,s,pt) = }_^nJ^2n{b,s) . (19) 

Tl = l 

More generally, for inclusive 2k jet production we have 



n>k 



J^2k{b,s)^Yl [Z]^2n{b,s). (20) 



This formula says that to every inclusive 2k jet produc- 
tion process there is a 2n jet exclusive contribution where 
we have combination of choosing k pairs from n pairs 
of jets. Equation (^01) can readily be inverted to find the 
exclusive Af2k{b, s) in terms of the inclusive J^2n{b, s), 

M2k{b,s) = U (-l)""'-^2„(6,s) . (21) 

n>k ^ 

The contribution r"gj](6, s) from minijet pairs to the to- 
tal inelastic cross section is by definition the sum over 
the individual exclusive contributions, 



r™t1(s,&) = E-^2fe(&,s). (22) 

k=l 

Using (|2ip in produces the simple expression, 

oc 

r"\1(s, b) = V(-l)"-W2„(6, s) . (23) 



n=l 



r*"j^(s,6) would be exactly equal to 



r"^'(s,6) = 2^(s,6)-|^(s,6)|^ 



if minijet production were present in all inelastic events. 
In general, therefore, we have the unitarity constraint, 

rjet^(s,o)<r (s,6). 

Of course, to evaluate r™fj(s, 6) one needs to know how 
to calculate the production cross section in impact pa- 
rameter space for arbitrary number of jets, i.e. one needs 
a model for each term in Eq. (PT|) . 



For this paper, we adopt the simplifying approxima- 
tions of Eqs. (Uni [H]) of simplifying (yet reasonable) as- 
sumptions. As already mentioned, we consider large im- 
pact parameters so that we may neglect correlations. If 
we interpret N'2{b,x,p^) as the probability to produce a 
minijet pair at impact parameter 6, then we may write 
the inclusive impact parameter cross section for produc- 
ing 2k minijet pairs as, 

M2k{s,b)^M2k{b.x,p'i) = {a\fMb,x,p1)f . (24) 

The expression (|23p for a given value of impact param- 
eter is then a geometric series that yields the analytic 
expression, 



r;"t1(s,6) = E(-i)""'-^^ = T^ 



a\f,^P 2{b,x,p'i) 
P2{b,x,pf, 



171C 

2'jet 



(25) 

If the final state partons are identical, then each term in 
Eq. should carry a symmetry factor of 1/fc!. In that 
case the series in Eq. ([25]) becomes. 



r™^(s,6) = l-exp[-a™:,J^(6,i,p^) 



(26) 



We will investigate both of these cases. Note that, as 
discussed in 9], Eqs. f [^5)1 and f [^5)1 satisfy the unitarity 
condition 



by construction. Equations ([25]) and (|26p are a reason- 
able approximation at values of b where N4{b, s) /N2{b, s) 
is a small parameter so that corrections due to correla- 
tions are small. It turns out that this is true for 6 ~ 2 fm 
for energies up to about ^/s « 50 TeV. 

In the specific case of four jets the profile function is 
defined as 



P,{b) 



P2{bf 



^(PbP2{bY 



where we have suppressed the other arguments in the 
profile functions. Note that P4 satisfies normalization 
condition by construction. Likewise the cross section for 
the four jet production is given by 



'2 jet 



(fbP2{bf 



B. Numerical Calculations 

Figure [3] shows sample calculations using Eq. (|25p for 
the case of (a): -^s = 14 TeV, relevant to the LHC, 
and for (b): -^/s = 50 TeV, relevant to cosmic ray air 
showers. Here we assume that the partons are identical 
particles. To test the sensitivity to the type of parton 
distribution being used, we have repreated the calcula- 
tion using MRST parton distributions in plots (c,d). We 
have done the calculation using two values of pt in each 
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FIG. 3: In all of these curves, the solid black curve is the extrapolation of the fit of Islam et al. (keeping only the diffractive, 
assymtotic limit). The black dotted curve is the calculation with pt = 2.5 GeV and the red-dashed curve is for pt = 3.5 GeV. 
The upper plots are calculated using CTEQ6M parton distributions and the lower plots are calculated used MRST parton 
distributions. (a,c) -/s = 14 TeV. (\),A)^/s = 50 TeV. These calculations are done assuming identical partons - Eq. lf26|l . 



case: the lower red-dashed curve is Wiih pt — 3.5 GeV, 
and the upper dotted curve is with pt — 2.5 GeV. In both 
of these calculations, we have used CTEQ6M gluon dis- 
tribution functions. For comparison, we have shown the 
profile for the total inelastic cross section obtained from 
the fits of Islam et al ^2^] (It is difficult at present to 
estimate the uncertainties in the extrapolation of mod- 
els for the inelastic impact parameter profile to the LHC 
energy. However the first data from LHC on elastic pp 
scattering will essentially eliminate this uncertainty.). In 
Fig. [3^, the calculation with a cutoff at pt = 2.5 GeV re- 
sults in a value for Vjlf^ that is larger than the inelastic 
profile function, F"*'^' at 6 > 1.5 fm, whereas the cut- 
off of = 3.5 GeV leads to a violation only at small 



values of b where the uncorrelated approximation likely 
breaks down. Therefore, the pt = 2.5 GeV cutoff is cer- 
tainly inappropriate for the maximum LHC energies. In 
Fig. we show the calculation for the higher energy 
y/s = 50 TeV. Here we see that the pt needs to be pushed 
even higher than 3.5 GeV to avoid a contradiction with 
unitarity. To test the sensitivity to the type of parton 
distribution being used, we have repreated the calcula- 
tion using MRST parton distributions in plots (c,d). We 
see that the above conclusions are essentially unchanged. 

It is worth emphasizing that typical values of b for 
jet production are much smaller, ^ 0.8 fm than the ones 
for which we obtained the constraint. Since the typical 
gluon densities involved for b 0.8 fm are significantly 
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larger (of the order of factor 3, compare Eq. (fT4|) ) than 
for 6 ~ 1.5 fm, the taming mechanism should be effective 
in the much broader range of pt ■ 

This result supports the cutoffs used in DPMjet II and 
III [3J in which the cutoff at this energy scale is roughly 
6 GeV, and SIBYLL 2.1 in which the cutoff is roughly 
5 GeV independent of b. However, in our logic it seems 
natural for the cutoff to inrease with a decrease of b. 

Our results disagree with the use of a low fixed cutoff as 
m QGSjet [i where the minimum pf is fixed at 1.5 GeV, 
and neXus 4] where the cutoff is fixed at 2.0 GeV. Note 
that it is to natural to expect that the pf cutoff is propor- 
tional to average gluon density encountered for a given 
impact parameter. Hence our limit for b ^ 1.5 fm cor- 
responds to a much larger value of the cutoff for more 
typical 6-values for the hard collisions - say production 
of two jets where average 6 are ^ 1 fm. 

We have also repeated the calculation using the as- 
sumption of non- identical partons - Eq. (|25p . This is 
shown in Fig. 01 We see that our above conclusions re- 
garding the minimum pt are not substantially affected. 

One source of uncertainty is in the true shape of the 
t-distribution of the two-gluon form factor. To test the 
sensitivity to this, we have repeated the analysis of the 
previous section using the more commonly used Gaussian 
model for the t-dependence of P2 rather than the two- 
gluon form factor. We fix the slope by demanding that 
(6^) for the Gaussian parameter equal the value of (5^) 
determined from the two-gluon form factor. The result 
is shown in Fig. [S] In this case, it seems that using the 
Gaussian parameterization leads to a greater violation of 
unitarity than the GPD. 

V. INTERACTION RADIUS AT VERY HIGH 
ENERGIES 

In this section we give a more general discussion of 
semi-hard interaction radius. We argue that it should be 
at least as large as what is used in the previous section to 
maintain consistency with DIS. A common approach to 
modeling the impact parameter dependence in extrapola- 
tions to very high energies is to write the profile function 
as, 

r(s,&) = l-cxp[-x(s,6)]. (27) 

Here x(s, h) is assumed to describe the phase shift pro- 
duced by each scattering of the proton constituents. 

In extrapolations to ultra-high energies it is also as- 
sumed that the basic parton-parton profile function can 
be decomposed into the sum of a term for soft scattering 
and for hard scattering: 

r(s,fe) = l-exp{-Xh(s,6)-Xs(s,6)} . (28) 

Corrections to this simple model are taken into account 
in more sophisticated versions of the eikonal model that 
include, for example, diffraction and triple Pomeron ex- 
change 0], though for the remaining discussions of this 



section, it will be sufficient to use Eq. ([28)1 . The soft 
part Xs{s, b) is modeled phenomenologically, whereas the 
semi-hard term Xh{s, b) can be calculated with the aid of 
Eq. However, we stress that since each term in the 
exponent of Eq. (|27|) is quite model dependent the rela- 
tionship between the unitarity of the total cross section 
and the large size of the semi-hard contribution, Xh{s, b) 
is unclear, and is likely to depend strongly on how we 
model the soft parts, etc. One can extract constraints on 
the UHE cross section indirectly by noting that the inclu- 
sive cross sections are obtained via sums over exclusive 
2?! minijet cross sections, 

oc 00 

Crtot = ^CT2n, <J2"iet — " ■ (29) 

71—0 n—1 

(Here it is assumed that the total cross section is domi- 
nated by jet pairs.) As in the previous section, we use a 
tilde on exclusive quantities. The total and elastic cross 
sections are, from Eqs. ^ and ([5]), 

atot =2jd^h[l- cxp{-xh(6) - Xsib)}] , (30) 

aeia = J d'h [1 - exp {-Xh{b) - xsmf . (31) 

In [l^, for example, a minimum value for p^ was deter- 
mined by requiring fits of (Ttot and ada. to data for the 
elastic cross section to be consistent with Eqs. ((29|) . In 
the semi-hard eikonal factor, a simple Gaussian model is 
used, 

X/.W = ^exp{-6V4i?g}. (32) 

Ro is the interaction radius in impact parameter space 
at the reference energy scale. (In general, there should 
also be a non-zero Regge slope a' leading to diffusion in 
the transverse plane as energy is increased.) 

A. Comparison with Deeply Inelastic Scattering 

From Eq. pil Bp and Eq. ([5^ we can extract the equiv- 
alent of Eq. (fTT)) for the Gaussian model, 

Pr"^^(5,.,pD = ^exp{-6V4i?g} . (33) 

When fitted to data the two over-lap functions Eq. ([55]) 
and Eq. (|17p should be in rough agreement. In particu- 
lar, both models should yield similar numerical values for 
(6^) calculated with Eq. For example, if Eq. §^ 

is used with a radius, Rq = 3.5 GeV~^ and a transverse 
momentum cutoff oi p^ = 3.5 GeV, one obtains a mean 
impact parameter of 0.75 fm, whereas if Eq. (fT7)) 

is used one obtains yJW) = 0.87 fm. In this case there 
is at least rough agreement between the two models. 
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FIG. 4: In all of these curves, the solid black curve is the extrapolation of the fit of Islam et al. (keeping only the diffractive, 
asymptotic limit). The black dotted curve is the calculation with pt = 2.5 GeV and the red-dashed curve is for pt — 3.5 GeV. 
The upper plots are calculated using CTEQ6M parton distributions and the lower plots are calculated used MRST parton 
distributions. Furthermore, the calculation is done assuming non- identical partons - Eq. (|25p . (a,c) y/s = 14 TeV. (b,d)v^ = 
50 TeV. 



In Ref. [10| it is found that a fit to the total cross sec- 
tion is also possible with a radius of Rq = 1.5 GeV~^ and 

= 2.5 GeV. However, in this case Eq. ([5^ produces 
a mean impact parameter equal to ( of 0.48 fm, 
whereas using Eq. (ini) produces ^/{W) = 0.89 fm. For 
this scenario, there is clearly an inconsistency between 
the Gaussian model and the impact parameter depen- 
dence extracted from DIS. The spread in impact param- 
eter space for these two scenarios is shown in Fig. [HI 

We can gain further insight by examining how the 
width of the impact parameter dependence of Tg{x^ p, /i) 
relates to the total qq-proton cross section. In DIS at 
high energy (or low-x) , cross sections are usually written 



in the target rest frame in term of the cross section for 
the interaction of a qq dipole with the target proton. The 
qq cross section is. 



"tot 



{d,x) 



-(fas{X/(f)xg{x,\/d?), (34) 



where d is the transverse size of the qq pair, and x is the 
longitudinal momentum fraction carried by the struck 
gluon 25| . Given the total cross section for the gg-proton 
interaction and the GPD in Eq. (fT5|) . we can write the 
amplitude for dipole scattering as, 



tid,x)Fg{x,t,fi) , 



(35) 
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4 and lOfTG*!. Because of the universality of the GPD, 
the gluon density in Eq. ([M]) should be the same as that 
appearing in Eq. ^ so long as we consider a dipole size 
d such that X/cP ~ p^. We are interested in the region 
of the integral close to and for typical values of x ~ x. 



Therefore, we evaluate Eq. ([55]) at dc = y^A/Jjjfj^. 
Using the FSW expression provided in Appendix B for 
mg{x,p'i) yields the solid curve shown in Fig.[71 The pro- 
file function calculated this way is well within the bounds 
of the unitarity constraint which is consistent with obser- 
vations at HERA. 



FIG. 6: A comparison of the distribution in impact parameter 
space for minijet production as predicted by the Gaussian 
model, Eq. (|33p and the 2-gIuon form factor, Eq. (|17|l . 

Since we are in the high energy limit, we neglect a small 
real part. is the center-of-mass energy of the dipole- 
proton coUins. We then invert Eq. [35]to obtain a profile 
function for the (/^-proton cross section for a pair of size 
d and for a gluon momentum fraction, x: 

rf(6,x) = 

aZ{d,^H{x,,) ^m,ix,,)b^^ K,im,ix,m. (36) 

In Eq. ([2]), the hard scale at which the gluon density 
and the strong coupling are evaluated is pt whereas in the 
DIS expression, Eq. it is related to the inverse size 
of the dipole's transverse size, X/d^, where A is between 



However, for (6^) to agree with the prediction from the 
Gaussian model with Rq = 1.5 GeV~^ andpf = 2.5 GeV, 
we find that we must increase the value of rrig by about 
a factor of 1.9. This leads to the dashed profile func- 
tion shown in Fig. [7| which exceeds the unitarity bound 
at small values of b. From Fig. [S] we see that, for the 
DIS calculation to remain consistent with unitarity, the 
Pt cutoff in the hadron-hadron jet cross section Eq. ([2]) 
should actually be pushed higher as the profile for hard 
partons becomes narrower. In other words, a narrow pro- 
file for hard scattering leads to a more rapid approach to 
saturation-like physics in the DIS calculation. To avoid a 
conflict with the unitarity in DIS, and to avoid a contra- 
diction with the observation at HERA that the amplitude 
for the dipole scattering is relatively far from the black 
disk limit even when x = 10"'' ~ 10""^ (at least unless the 
scales are greater than ~ 2 GeV^, see [l^) we must use 
a radius-squared for the hard interaction that is greater 
than - 3.5 GeV-^. 
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FIG. 7: The profile function (solid curve) calculated using 
Eq. 1361 To match with the Gaussian model, mg must be 
increased by a factor of 1.9 - resulting in the dashed curve. 



VI. CONCLUSIONS 



model which includes the correlations for multi-jet pro- 
duction in the impact parameter space. 

We have also analyzed in detail the dominant regions of 
the integrand in the minijet formula and shown that the 
dominant configurations for small energies are when the 
jets are produced with approximately equal but opposite 
rapidities. When the energy is increased the configura- 
tions with the jets going into the same rapidity hemi- 
sphere are equally important. This includes the forward 
regime in hadron scattering where the saturation correc- 
tions in the gluon density are likely to be important at 
the LHC energy. 
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In this paper we have analyzed the region of appli- 
cability of the pQCD formula for the minijet production 
in hadron-hadron collisions. Based on the unitarity argu- 
ments in impact parameter space we have set a constraint 
on the minimum value of pt for which the formula can 
be used. The model satisfies the unitarity by construc- 
tion, i.e. by the summation of the multiple scatterings, 
without any correlations in impact parameter space. For 
the LHC-scale energy ^/s = 14 TeV the minimum value 
of > 2.5GeV and for the cosmic ray energies for pp 
collisions = 50 TeV, pt > 3.5 GeV. 

We again stress that these values are for large impact 
parameters ~ 1.5 fm. Though our simple analysis is not 
effective at small impact parameters due to the neglect 
of correlations and due to r(s, b) ^ 1, it is reasonable to 
expect the minimum pi to rise as the gluon density rises. 
Thus, the true minimum pi is probably even larger than 
what we find here. 

For the case of the cosmic ray interaction of protons 
with air, the typical gluon densities are at least factor 
five larger than the gluon densities encountered in the pp 
case at 5 1.5 fm. 

Note that the values of the minimum pt are much larger 
than any conceivable soft scale. Therefore, the dynamical 
mechanism for violation of the leading twist formalism 
which we find here must include some kind of non-linear 
strong field effects at relatively small coupling. 

The analysis can be extended for a more sophisticated 



Appendix A 

Parametrization of the mass parameter rrig in the 
profile function 

The mass parameter nig in p3p is related to the inverse 
of the average impact parameter corresponding to the 
profile distribution (fT4|) 




In [2TI the following parametrization was found for the x 
and scale dependence of irig 

{p^){x,Ql) = max|o.31fm^-K0.0194fm2ln^,0.28fm^| 
and 

{p'){x,Q') = {p')ix,Ql)(^l + A\n^^ \ 

where 

Qg = 3GeV^ ^ = 1.5, a = 0.0090 In i. 
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